Lattice normal modes and electronic properties of the correlated metal LaNi0 3 
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We use density functional theory (DFT) calculations to study the lattice vibrations and electronic 
properties of the correlated metal LaNiOs • To characterize the rhombohedral to cubic structural 
phase transition of perovskite LaNiOs, we examine the evolution of the Raman-active phonon modes 
with temperature. We find that the A\ g Raman mode, whose frequency is sensitive to the electronic 
band structure, is a useful signature to characterize the octahedral rotations in rhombohedral LaNiOs- 
We also study the importance of electron-electron correlation effects on the electronic structure with 
two approaches which go beyond the conventional band theory (local spin density approximation): 
the local spin density+Hubbard U method (LSDA+?7) and hybrid exchange-correlation density 
functionals which include portions of exact Fock-exchange. We find the conventional LSDA accurately 
reproduces the delocalized nature of the valence states in LaNiOs and gives the best agreement to 
the available experimental data on the electronic structure of LaNiOs. Based on our calculations, we 
show that the electronic screening effect from the delocalized Ni 3d and 0-2p states mitigate the 
electronic correlations of the d 7 Ni cations, making LaNiOs a weakly correlated metal. 

PACS numbers: 71.15.Dx,63.20.dk,78.30.Er,79.60.Bm 



I. INTRODUCTION 

Conducting electrode materials are critical elements in 
the design of ultra-thin ferroelectric devices magne- 
toresistive elements P and magnetoelectric multiferroic 
memories^ The performance of perovskite-based het- 
erostructures is intimately related to the electronic and 
atomic coherency across the electrode-oxide interface. 7 
For this reason, suitable metallic perovskite oxides that 
are structurally compatible with the active functional 
layers are highly desirable. LaNiOs (LNO) has found 
widespread use in this capacity because it is a 3d 
transition metal oxide that shows no metal-insulator or 
spin state transitions! 1 1 * 12 l At present, it is also actively be- 
ing pursued as the functional oxide component in devices 
that could harness an electric-field tunable Mott insula- 
tor transition (so called "Mottronic" applications)}^^ 
because it is the end-member of the rare-earth nicke- 
late series — compounds with small charge-transfer gaps 
that are highly susceptible to temperature^ and pressure- 
induce d 17 ! 18 ! electronic phase transitions. 

Metallicity and magnetism in 3d transition metal (TM) 
oxides are strongly dependent on the valence bandwidth, 
which originates from the hybridization between the TM 
3d and O 2p orbitalsP^In perovskite oxides like LNO, the 
hybridization derives from the structural connectivity of 
the Ni06 octahedral units throughout the crystal. Small 
changes in the Ni-O-Ni bond angles and the Ni-0 bond 
lengths, therefore, can dramatically alter the electronic 
properties. For example, a tunable insulator-metal (IM) 
transition is obtained from the isovalent substitution^ 
of La with rare-earth elements: The charge-transfer gap 
closes and the IM-transition becomes accessible above 



room temperature as the rare-earth ionic radius increases 
and straightens the Ni-O-Ni bond angle. 

In addition to changes in the structural degrees of 
freedom, the small spatial extent of the 3d orbitals also 
reduce the valence bandwidth. This sufficiently enhances 
electron-electron correlation effects such that conventional 
band metals are often rendered insulating In LaNiOs, 
however, the strong Ni 3d - O 2p covalent interactions — 
formally it contains Ni 3+ cations in a low-spin 3d 7 con- 
figuration (t2 q e )j) — are anticipated to reduce the corre- 
lation effects) 2 ^ Nonetheless, there are clear signatures 
indicative of important electron-electron interactions in 
the T 2 depen dence of the resistivity and heat capacity 
measurement sPEU Magnetic susceptibility data also re- 
veal en hance d Pauli paramagnetism and effective carrier 
masses! 16 * 23 " 26 ^ Consistent with those studies, temperature 
dependent x-ray photoemission finds that the spectral 
weight of the Ni e g band at the Fermi level increases upon 
cooling! 2 ^ A fundamental question regarding the intrinsic 
properties of LaNiOs therefore still remains: which fac- 
tors of the interwoven (correlated) electronic and atomic 
structure support its metallic state? 

In this work, we perform first-principles density- 
functional theory (DFT) calculations to investigate how 
hybridization between the Ni 3d and O 2p states and 
the structurally correlated Ni06 octahedral framework re- 
spond to electron-electron interactions. We first examine 
the temperature-dependent rhombohedral to cubic struc- 
tural phase transition using a Landau formalism 28 and 
ab initio- derived phenomenological coefficients obtained 
from the conventional band theory [local spin density 
approximation, (LSDA)]. In the rhombohedral phase, we 
compute the evolution of the Raman-active phonon modes 
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with temperature and find the LSDA results give good 
agreement with experiment. We then show how the A\ g 
Raman mode, which describes the rotation of adjacent 
Ni06 octahedra, can be used as a structural indicator for 
the low-temperature rhombohedral phase. 

To study the coupling between the lattice modes and 
the electronic structure, we compare conventional LSDA 
results with three other approaches designed to improve 
accuracy: the local spin density +Hubbard U (LSDA+C/) 
method, and t wo hybri d exch ange- correlation function- 
al, PBECP^Hland HSEP^ which contain mixtures of 
Fock-exchange added to the generalized gradient corrected 
DFT-functional of Perdew, Burke and Ernzerhof (PBE)P 
We then examine the electronic structure of LaNiC>3 by 
comparing the results obtained from the various func- 
tional with the experimental photoemission spectroscopy 
(PES) and x-ray photoelectron spectroscopy (XPS) data. 
We find that the screening effects originating from the 
hybridized O 2p and Ni 3d electrons are sufficiently strong 
that they reduce the electronic correlations in LaNiOs, 
making it a weakly correlated metal. 



II. CRYSTAL STRUCTURE AND RAMAN 
MODES 

The crystallographic tolerance factor 3 - 4 ^ for LaNiOs is 
£=0.97 and indicates that the aristotype cubic phase is 
susceptible to octahedral rotations because of an under- 
bonded La-0 coordination.^ At ambient conditions, bulk 
LaNi03 crystallizes in a rhombohedral structure with 
space group R3c [Fig. [ija)] and is related to the cubic 
perovskite through a trigonal lattice distortion along the 
[111] -body diagonal that doubles the primitive unit cellP^ 
It also exhibits octahedral rotations, which are equal in 
magnitude and alternate in "sense" about each Cartesian 
direction; this tilt system is classified according to Glazer 
notatio n 1 35 * 36 ^ as a~a~a~ [Fig. [I] (b)]. Similar to other 
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FIG. 1: (Color) (a) The crystal structure of rhombohedral 
R3c LaNiOs possesses anti-phase rotations (a~a~a~ tilt sys- 
tem) of adjacent NiC>6 octahedra with angle (b) about the 
[lll]-trigonal axis. The relationship between the pseudo-cubic 
perovskite cell (solid lines) and the rhombohedral lattice vec- 
tors (dashed lines) of length a is also illustrated. 
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FIG. 2: (Color) Illustration of the displacement patterns for 
the Raman-active lattice normal modes in the R3c structure. 
The corresponding symmetry labels and our calculated LSDA 
frequencies (in cm -1 ) are given for reference. 

rhombohedrally distorted perovskitesp^ LaNiOs under- 
goes a temperature-induced rhombohedral-to-cubic phase 
transition 38 upon heating: the Ni06 octahedral rotation 
angle reduces continuously to zero in the cubic phase. 

In the R3c space group, the La cations occupy the 
2a (|, |, |) Wyckoff positions while the Ni cations occupy 
the 2b (0, 0, 0) positions. The oxygen atoms are at the 
6e (x,x+ \ , \) site, where x is the only free internal struc- 
tural parameter that sets the rotation angle of the Ni06 
octahedra. These 10 atoms in the primitive rhombohedral 
unit cell give rise to 30 zone-center vibrational modes with 
the irreducible representations (irreps): 
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• 4E g + 2A lu + 4A 2u + 6E U . 



where the E g and E u modes are two-dimensional irreps. 
The infrared-active and acoustic modes transform as the 
4^2n + 6£? n , while the Raman active modes are given as 
Ai g + 4:E g . The A\ g and E g Raman modes of the rhom- 
bohedral structure are mainly related to collective modes 
of the oxygen octahedral network (Figure |2j: The A\ g 
mode describes the rotations of the Ni06 about the trig- 
onal [lll]-axis, and the E g modes describe anti-parallel 
La displacements, Ni-0 bond bending, stretching, and 
octahedral rotations about axes perpendicular to the [1 In- 
direction. 



III. THEORETICAL METHODS 
A. LSDA 

We use two DFT implementations in this study, the 
quantum ESPRESSO code (QE] P an d the Vienna 
ab initio Simulation Package (vasp)! 40 ! 41 * Our reference 
electronic structure model to which we compare our ad- 
vanced treatments of correlation effects is the frequently 
used local-spin density approximation (LSDA) exchange- 
correlation (XC) functional! 4 ^ Both DFT codes use the 
Perdew-Zunger parametrization 42 of the Ceperley- Alder 
datsF 3 ^ for the XC-functional. In all calculations, we con- 
strain a collinear spin configuration for the Ni ions. The 
atomic positions in the rhombohedral structure are opti- 
mized by starting from the positions reported in Ref fTT| 
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and the ionic coordinates are relaxed until the Hellmann- 
Feynman forces on the atoms are less than 0.1 meV A -1 . 

VASP Details. We treat the core and valence electrons 
for all calculations in VASP using the projector augmented 
wave (PAW) methocF 4 -' with the following valence elec- 
tron configurations: bp 6 bd 1 6s 2 (La), 3p 6 3d 9 As 1 (Ni), and 
2s 2 2p 4 (O) and a 650 eV plane wave cutoff. We use a Gaus- 
sian smearing of 0.05 eV over a 7 x 7 x 7 Monkhorst-Pack 
/c-point mesh 45 centered at T for the Brillouin zone (BZ) 
integrations (172 points are sampled in the irreducible 
BZ) in the 10-atom rhombohedral unit cell. 

We obtain the phonon frequencies at the T-point by 
calculating total energies with respect to atomic displace- 
ments from the reference R3c structure. In this frozen- 
phonon method a series of small (symmetry inequivalent) 
atomic displacements are imposed along different Carte- 
sian directions. We calculate the dynamical matrix from 
the Hellmann-Feynman forces 4 ^ induced on the ions after 
making the small positive and negative displacements 
(to remove any quadratic effects) about the reference 
structure positions from total energy DFT calculation. 
Diagonalization of the dynamical matrix yields the atomic 
displacement patterns (eigenvectors) and phonon mode 
frequencies (eigenvalues). 

QE Details. In our plane-wave calculations with the 
QE code, we use norm-conserving^ optimized designed 
nonlocaP^ pseudopotentials. The following valence elec- 
tron configurations are used: 5s 2 5p 6 4/°5d 6s 6p°(La), 
3d 9 4s 1 Ap° (Ni) and 2s 2 2p 4 (O). Note, In both VASP and 
QE, partial core corrections (PCC) are included in the 
Ni pseudopotentials following the prescription in Refl49l 

We sample the BZ using an 8 x 8 x 8 Monkhorst-Pack 
fc-point mesh. A plane-wave energy cutoff of 50 Ry was 
used for calculation. Good convergence can be reached, 
as test calculations with higher energy cutoff (up to 75 
Ry) yield the same results. We obtain phonon modes and 
displacement vecto rs usin g density functional perturba- 
tion theory (DFPT)) 50 * 51 * by constructing the interatomic 
force constants from the dynamical matrices obtained on 
a uniform 4x4x4 g-point grid. The phonon frequencies 
and vibrational eigenvectors at arbitrary q vectors are 
then calculated from diagonalization of the dynamical 
matrix. 

We point out here that both implementations of DFT 
at the LSDA level give good structural agreement (Table 
[l|. This fact is important for the comparison study we 
make later which uses various functionals available in the 
different codes. Each code yields a small underestimation 
of the experimental lattice constant which is typical for 
the LSDA. The QE code more closely reproduces the ex- 
perimentally measured Ni-O-Ni rotation angles, whereas 
the VASP code gives a better rhombohedral angle (a r ho)- 
For the electronic and magnetic structure of LNO, both 
codes predict LNO as a non-magnetic metal, which is the 
reasonable representation of the experimentally observed 
Pauli paramagnetic ground- stated Detailed electronic 



rect ground-state of paramagnetic LNO. Without PCC, 
LSDA calculations yield a spin-polarized solution with a 
non-zero local magnetic moment on the Ni cations. 



B. LSDA+J7 

A complete description of correlated 3d transition metal 
oxides with narrow valence bandwidths is challenging 
within a density functional approach due to the limita- 
tions of standard exchange-correlation potentials in de- 
scribing localized electronic states. To remedy this prob- 
lem, we use "beyond-LSDA" techniques, starting with the 
LSDA+Z7 approach. This method is generally regarded to 
be the most computationally feasible means to reproduce 
the correct ground-states in correlated systems.^ In this 
formalism, the LSDA energy functional is expanded to 
include an additional on-site orbital-dependent energy 
term cast as a Hubbard repulsion U and an intratomic 
Hund's exchange energy J. 

To reduce the ambiguity in the definition of the 
LSDA+/7 parameters, we use the spherically averaged 
form of the rotationally invariant LSDA plus Hubbard 
U method^ introduced by Dudarev et al. with only 
one effective Hubbard term, f7 e ff = U — J. We treat the 
double-counting term within the fully localized limit. 54 

The change in total energy Ejj for including the Hub- 
bard correction to the exchange-correlation potential is 



r), (1) 



where n ma are the spin (a) and orbital (m) occupation 
numbers at site i. In the limit where the occupation 
matrices are integer and diagonal, the LSDA+/7 correc- 
tion can be understood as a shift in the occupied orbitals 
(n m = 1) by — U/2 to lower energy and by -{-U/2 higher 
in energy for unoccupied orbitals (n m = 0). 

In this study we examine U e s values of 0, 3 and 6 eV for 
the correlated Ni 3d orbitals states; the standard LSDA 
corresponds to a Z7 e ff = eV. We note that throughout 



TABLE I: Structural parameters obtained for R3c LaNiOs 
within the LSDA from the two codes used in this work and 
their comparison with experimental data (cf. Ref|Tl ). In 
the rhombohedral setting, the rotation angle of adjacent Ni06 
octahedra along the trigonal axis is given as = arctan(2y / 3'^), 
where x = u + § . 
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properties for LNO will be presented in Sec. IV Here 
we note that PCC for Ni is crucial to obtain the cor- 





Experiment 


VASP 


QE 


X 
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FIG. 3: (Color) The occupation numbers n* for the total Ni 3d, 
e g and tig orbitals in LaNiOs, with respect to potential shift a 
for the initial (0) (filled circles) and self-consistent (SCF) (open 
squares) electronic configurations. The slopes of the linear 
data fits are used to construct the %o and x response matrices. 
The occupation numbers are normalized by substracting off 
the occupation values m(0) with a zero potential shift (a = 0). 



the remainder of the manuscript, U denotes the effec- 
tive Hubbard parameter. We will discuss the structural 
changes induced by varying U in the subsequent sections. 



1. Self- consistent Hubbard U 

We also calculate a self-consistent U parameter for 
LaNiC>3 following the scheme developed by Cococcioni 
and de Gironcoli. 55 Their approach relies on the fact that 
the potential U restores the correct piece-wise linear be- 
havior of the system's total energy as function of electron 
number, whereas the LSDA functional incorrectly predicts 
parabolic dependence of energy on occupation number. 
The effective interaction parameter U is deduced from 
Janak's theorem and linear response theory as 



U = (Xo 1 -X ^iijwithxi. 



dn l 
da. 



(2) 



Here n l is the occupation number of the localized levels 
at site i and otj represents the potential shift applied 
on the localized orbital at site j. The response matrix 
Xo describes the noninteracting contribution to the band 
structure after application of a potential shift a, while x 
represents the fully self-consistent response to the same 
potential shift. In practice, the first term xo is com- 
puted from the first iteration in the self-consistent (SCF) 
electronic minimization. 

We compute the linear response of the occupation num- 
ber n l using the LSDA functional and norm-conserving 
pseudopotentials within the QE package. The full re- 
sponse matrices xo and x are then computed by perform- 
ing the linear response calculation within a 2 x 2 x 2 
pseudo-cubic supercell (40-atom cell, containing 8 Ni 
atoms), which is large enough to give a converged U value 
and avoid spurious interactions from the local potential 
a on neighboring Ni sites. 

The low-spin Ni 3d orbital occupation anisotropy — a 
fully filled ti g manifold and a quarter-filled e g shell — 



suggests that the effective Coulomb repulsion experienced 
by an electron on different orbitals should also be substan- 
tially different For this reason, we further decompose 
the orbital occupancy n l — » (m is the angular mo- 
mentum quantum number, and orbital occupation 
here refers to Lowdin atomic charge) into crystal-field 
split (i-manifold contributions by projecting the valence 
states on to atom-centered sites i with a Lowdin 

orthogonal atomic basis set as 

vu\n v ) = \ EE (! - ^<a)Wm)mn v ) ■ (3) 

i ma 

By replacing the orbital occupancy n % in Eq. (2) with 
n l m , we can calculate the effective Hubbard term U for 
the two independent Ni-c? manifolds: U(e g ) and Ufog) 
respectively. Since the t<i Q manifold is fully occupied 
\ n (p2g) ~ 1], the linear response calculation of Ufag) 
requires that the contribution from these states to the 
total Coulomb interaction be nearly zero. Within the 
numerical noise of our calculations we find this to be the 
case. Fig. 3 shows that difference in the xo and x response 
functions (slopes) of the t^g states is practically negligible. 
The Ni e g states, however, which are partially filled have 
a substantial non-zero interaction, U(e g ). These states 
lead to a self-consistent Hubbard U value of 5.74 eV for 
LaNi0 3 . 

Recently Nohara et at. studied LaNiOa with LSDA+Z7 
and GW methods and fitted their calculated energy 
spectra to experimental XPS and x-ray adsorption spec- 
troscopy (XAS) datsF^ to obtain a Coulomb interaction 
of U = 7.0 eV and an exchange interaction J = 1.3 eV for 
the Ni d-manifold. These values correspond to an effective 
Hubbard U term (f7 e ff = 5.7 eV) which is consistent with 
our linear-response calculation of U. 



C. Hybrid Functionals 



Recent studied ^ 33 * 59 1 have demonstrated that local or 
gradient-corrected density functionals for exchange do 
not closely reproduce Hartree-Fock calculations, and that 
inclusion of some exact Fock (F) exchange improves elec- 
tronic structure properties including band gaps, orbital 
localization and electronic polarizations.^ Accordingly, 
in the present work, we use the PBE0 and HSE hybrid 
functionals as implemented in the QE package and then 
contrast them to the LSDA. 

The exchange terms (V x ) in the PBE0 functional are 
constructed by mixing 25% of exact-exchange with 75% 
GGA-PBE exchange,^ while the electron correlation (V c ) 
part is represented using only the correlation components 
from the PBE functional 33 : 



PBE0 



1 



3 yPBE 

4 x 



PBE 



(4) 



The one- quarter mixing parameter of exact exchange 
in the PBE0 functional was obtained by fitting to the 
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atomization energies of a large database of molecules.^ 
For periodic bulk systems, however, it is argued that the 
best fraction of exact exchange for modeling s olid-st ate 
electronic structure is highly system-dependent ! 29 l 61 l 

The HSE hybrid functiona l 31 * 32 * is less computationally 
demanding than PBEO (unless the recent Wannier func- 
tion implementation described in Ref.30 is used), since 
it avoids the slowly decaying Fock-exchange interactions 
by substituting part of the long-range Coulomb-kernel 
with a density functional analog. This approximation 
makes the HSE functional more efficient for plane-wave 
calculations of periodic systems. The expression for the 
exchange-correlation potential is given as 

yHSE = 1 ysr,» + 3 y PBE,sr, M + yPBE.lr.M + y PBE > (5) 

where \i is the parameter that controls the decomposition 
of the Coulomb kernel into short-range (sr) and long- 
range (lr) exchange contributions. This type of screened 
exact-exchange is absent in the PBEO functional. In this 
work, we follow the HSE03 parametrization 31 which sets 
the cut-off distance to \i = 0.109 A -1 . 

Similar to PBEO, the Fock-mixing parameter in HSE 
varies approximately in proportion to the inverse static di- 
electric constant (s^ 1 ), and therefore also becomes highly 
system dependent. For metals with excellent screening, 
approaches zero. Therefore PBEO or HSE calculations 
which incorporate some Fock-exchange may lead to severe 
overestimation in orbital bandwidths and spin-exchange 
splitting parameters for such materials.^ 

In our hybrid functional calculations on metallic 
LaNiOs, we use a dense Brillouin zone sampling (8x8x8 
/c-point grid) for accurate evaluation of the Fock-exchange 
operator in reciprocal spaced Due to the heavy compu- 
tational cost, we are unable to perform the structural 
optimization for LNO using the hybrid functionals in the 
present work. According to our test calculations using 
the LSDA, PBE and LSDA+Z7 functionals, the subtle 
differences between the ground-state atomic positions 
produce negligible changes to the computed electronic 
band structure. We therefore choose the LSDA optimized 
atomic structure as the ground-state structure for use in 
the hybrid functional calculations. 

IV. RESULTS AND DISCUSSION 

In this section, we present our first-principles results 
on the structural, vibrational and electronic properties of 
bulk LaNiOs obtained using the various XC-functionals. 
This section is organized as follows: we first investigate the 
structural phase transition and temperature-dependent 
Raman phonon modes of LaNiOs using the conventional 
LSDA approximation. We then examine how the atomic 
structure, Raman phonon modes and electronic prop- 
erties of LaNiOs are modified with various treatments 
of electron-electron correlation effects. Finally, by com- 
paring the densities-of-states results obtained from the 



different functionals to the experimental spectroscopic 
data, we identify which XC-functional best reproduces 
the intrinsic electronic properties of LaNi03. 

A. Structural Phase Transition 

The temperature-induced rhombohedral-to-cubic phase 
transition in LaNi03 manifests as a cooperative rotation 
of Ni06 octahedra along the trigonal lattice axis.^In this 
section, we show how the structural transition to the rhom- 
bohedrally distorted perovskite phase is characterized by a 
specific soft mode which has a frequency that collapses to 
zero as the system approaches the cubic phased We note 
that all results in this section are obtained with the LSDA 
XC-functional and NC-pseudopotentials as implemented 
in QE. 



1. Minimal Landau Model 

In the Landau theory framework for phase transitions, 
the free energy of the system is expanded in powers of an 
order parameter that characterizes the transition. From 
our previous discussion of the low-temperature R3c crys- 
tal structure of LaNi03, the only free internal parameter 
is the oxygen position x at the 6e Wyckoff site. There- 
fore, the natural order parameter to characterize the 
rhombohedral-to-cubic phase transition is a structure- 
adapted form of the free Wyckoff position: the NiC>6 
octahedra rotation angle (Fig. [T]). 

We expand the free energy G in even powers (up to 
fourth order) of this rotation angle order parameter as: 

G{0, T) = G (T) + k(T - T c )0 2 + X0 4 (6) 

where n and A are temperature-independent coefficients 
and is the angle of rotation about the [11 Indirection. 
We next "freeze-in" the NiC>6 rotation pattern [inset of 
Fig. 5 (b)] which corresponds to the A\ g Raman- mode of 
the R3c phase. In Fig. [3J we show our calculated values of 
the total energy versus Ni06 rotation angle computed at 
the LSDA equilibrium volume. Our fit of the data to Eq.[6] 
yields excellent agreement and approximately corresponds 
to the free energy at zero temperature G(0,T = K). 
From the minimum of the free energy (dG/dO)T — 0, we 
obtain the T dependence of 9: 

e 2 = ^(T G -T),iorT<Tc, (7) 

and a critical rotation angle at T = K of 0c = 10.08° 
(the optimized 6 from the LSDA calculation). 

We next calculate the energy stabilization obtained 
from freezing-in the NiOe rotation at K as the differ- 
ence between cubic and rhombohedral phases, AE = 
177 meV/10-atom unit cell given by the well depth from 
the total energy calculations [Fig. Qa)]. With increasing 
temperature, the energy stabilization from the structural 
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Ni0 6 rotation angle 6 (degrees) temperature (K) 

FIG. 4: (Color) (a) Landau free energy G(0, T = K) of 
LaNiOs as a function of the order parameter 0. Solid lines 
are calculated using Landau theory and the filled symbols 
correspond to the LSDA total energy results, (b) Equilibrium 
order parameter as a function of temperature; the second 
order phase transition occurs at T — Tc- 



distortion decreases until at T = Tc, the thermal exci- 
tation energy is equivalent to AE and the double-well 
potential becomes a single well with one minimum at 
= 0°. With these conditions, we write 

AE = (T • AS) T =T C , (8) 

where AS refers to the entropy difference between the cu- 
bic and rhombohedral phases. Since LaNiOs has a Debye 
temperature 0£> of 420 Kp^at the structural transition 
T = Tc, Tc ^> 0£>, we are able to confidently treat 
each lattice mode as an independent harmonic oscillators 
which subsequently contributes NksT to the free energy. 

For the 10-atom rhombohedral LaNiOs unit cell, a sin- 
gle soft mode (N = 1) describes the transition into the 
high-temperature cubic phase; therefore, AS — From 
these conditions, we calculate the LSDA critical tempera- 
ture for the structural phase transition Tc to be 2057 K, 
which is close to experimental result (Xc=1780 K). The 
overestimation of Tc likely originates from the calculated 
enhancement of #c, which increases the energy difference 
between cubic and rhombohedral phases. 

Using this calculated Tc value, and combining with 
our earlier calculated DFT total energy results, we obtain 
the coefficients in Eq. (6): «=1.696 /ieV/K-(°) 2 and A = 
17.1 /ieV/(°) 4 . These values combined with Eq. [7] give 
the second-order temperature dependence of the Ni06 
octahedral rotational angle 0(T) shown in Fig. Qb). 

2. Correlation of volume expansion with phase-transition 

With increasing temperature, the LaNiOs lattice under- 
goes a thermal volumetric expansion process. To precisely 
simulate the temperature induced rhombohedral-to-cubic 
phase transition process, it is necessary to evaluate the ef- 
fect of volume expansion on the Ni06 octahedral rotations 
across the structural phase transition. 

By including the coupling between the order parame- 
ter and the equilibrium atomic volume V, we extend 



the minimal ID Landau model discussed in the previous 
subsection into a 2D case, with the free energy G given 
as: 

G{9, V, T) = G (V, T) + k(V)(T - T c )0 2 + \{V)0 A , (9) 

where V is the equilibrium LaNiOs volume at a given 
temperature, and the coefficients k and A are expressed 
as a function of V. From K to room temperature, we 
assume a linear volumetric thermal expansion: V(T) = 
Vo x (1 + ayT), where Vo is the volume at K and the 
thermal volumetric expansion coefficient ay = 1.624 x 
10 _5 /K is from experiment.^ 

To correlate the DFT total energy results with the 
Landau free energy at K, we rewrite Eq. [9] as 

G(<9, V, T = K) = G (TO + A(V)0 2 + B(V)0 4 . (10) 

This form indicates a practical route to calculate the 
optimal NiOe rotation angle 0c with respect to V. We 
apply the procedure described in the previous subsection, 
but now obtain volume-dependent coefficients A and B by 
fitting each coefficient through a Taylor expansion about 
the equilibrium volume: 

AV AV o 

A(V) = A + A 1 (-f)+A 2 (-f ) 2 + ..., (11) 

AV AV o 

B(V) = B + B 1 (-^- ) + B 2 (-f- ) 2 + ... 

where AV = V — Vo- With this method, we obtain the 
following coefficients: Aq = —3.6, A\ = —16.6, A^ = 
1600.1 meV/(0 2 , and = 17.9, B x = -15.5, B 2 = 
-1664.0 /ieV/(0 4 . 

This treatment of changes in the lattice volume at the 
rhombohedral-to-cubic transition reveals that the magni- 
tude of our phenomenological coefficients A (B) increase 
(decrease) with increasing cell volume. Thus, the curva- 
ture of the double- well potential [Fig. [I] becomes steeper 
and its depth deeper as the cell volume increases; this fact 
indicates that the thermal expansion effect will hinder 
the rhombohedral-to-cubic phase transition in LaNiOs. 
Nonetheless, the extent of the thermal volumetric ex- 
pansion effect is very limited, i.e. from K to room 
temperature, AV/Vq ~ 0.5%, which leads to changes in 
the coefficients A and B by less than 1% and a devia- 
tion in our simulated rotation angle at 300 K by 0.2 ° 
compared to the volume-independent model previously 
described. We conclude that, for the temperature range 
we investigated (0 K to room temperature), treating the 
phenomenological coefficient in our Landau theory as 
volume-independent suffices to produce an accurate de- 
scription of the temperature-induced octahedral phase 
transition in LaNiOs . 

3. Temperature- dependent Raman Frequencies 

We next study the structural and vibrational properties 
of LaNiOs at finite temperatures. As the temperature 
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FIG. 5: (Color) (a) Temperature-dependent Raman frequen- 
cies for the rhombohedral LaNiOs structures. Solid lines 
indicate experimental results which are taken from Ref 63] Re- 
sults from our LSDA calculations are given as the filled circles, 
(b) We observe linear scaling behavior between Ai^ Raman 
frequency and the NiOe octahedral rotation order parameter 
0. The fitted solid line is a guide to the eye. The inset (upper 
left) shows the vibrational pattern of Ai^ mode (both side 
and top view), while the other inset (lower right) shows a 
linear change in the squared Ai g soft-mode frequency with 
temperature. 



increases the Ni06 rotation angle decreases [Fig. Qb)] 
and the Ni-0 bond length increases. To capture this 
effect in our simulation, we choose to combine the cell 
parameters a and a for LaNiOs at a given temperature 
with a specific NiC>6 rotation angle and Ni-0 bond 
length [d(Ni-O)]. The rotation angle for a given tem- 
perature is obtained from Eq. [7| While for <i(Ni-0), we 
refer to the experimental temperature-dependent neutron 
scattering results.^ To sample the d( Ni-O) -temperature 
space at experimental values not available, we assume 
linear thermal expansion. We then relax only the internal 
coordinates and use the structures that yield the best 
experimental agreement to represent a snapshot of the 
experimental LNO structure. 

We use these structures to perform phonon calculations 
and obtain the temperature-dependent Raman frequen- 
cies [Fig. J5^a)]. A comparison of the experimental data 
to our DFT calculated temperature-dependent Raman 



frequencies reveals good agreement: both data sets show 
red-shifts in the Raman frequencies with increasing tem- 
perature. The calculated high-frequency E g mode, which 
corresponds to stretching of the Ni-0 bonds, however 
is systematically underestimated by around 10%. The 
LSDA also predicts that the frequencies of the bending 
and stretching E g modes decrease nearly twice as fast as 
that experimentally observed. 

We now connect the frequency of the A i g Raman mode 
(<UA lg ) and the order parameter described in the previous 
subsections. Here, we assign the lattice mode to the 
stiffness k, or the curvature (second-order derivative) of 
the potential well, as 



U M 9 



oc 



d 2 G 
80 2 



4k(T c - T), for T<T C . 



(12) 



12 



we find that uja 1 oc 
T T c both and 



By comparing Eq. [7] with Eq. 
|0| and in the temperature limit 
WA lg approach zero. Using our fitted value of we 
show in Figure f5^b) that 0JA lg varies linearly with 0. 
We also predict from the slope of Fig. |5^b) that the 
Raman-mode shift due to changes in the NiOe rotation 
angle is uoa 19 /0 = 23.0 cm _1 /(°). Finally, we plot the 
temperature-dependence of the squared A\ g frequency 
[Fig. [HJinset)] and find that ^\ lg decreases linearly with 
temperature, confirming the uj\ {T) relation in Eq. (12). 
On the basis of these findings, we conclude that the A\ g 
soft mode is an excellent signature for the magnitude of 
the octahedral rotations in rhombohedral LaNiOs and its 
deviation from cubic phase. 



B. Correlation effects on the atomic structure and 
Raman phonon modes 

We focus in this section on how the Raman active 
mode frequencies are modified through changes in electron- 
electron correlations. We first decompose the effect of 
correlation through the Hubbard U term on the structural 
degrees of freedom by fixing the lattice parameter to 
that of the experimental R3c structure and allowing the 
internal atomic positions to fully relax. The results of our 
atomic relaxations for LSDA+/7 values of 0, 3, and 6 eV 
are summarized in Table \H\ In all cases, the LSDA+L 7 " 
functional accurately reproduces the known experimental 
ground-state structure with the minor caveats we discuss 
next. 

With increasing correlation, we find that the Ni06 rota- 
tion angle increases beginning from the LSDA structural 
ground- state (U = eV), which slightly underestimates 
the rotation angle, to U = 3 eV which overestimates it 
by approximately 1°. By further increasing the Hubbard 
U value to 6 eV, we find the rotation angle decreases. 
The consequence of keeping the unit cell volume and 
rhombohedral angle fixed is that the change in the NiOe 
rotation angle must be accommodated by bond stretch- 
ing (or compression) rather than through rigid rotations 
(constant Ni-0 bond lengths). Because the Ni atoms also 
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FIG. 6: (Color) Spin- and atom-resolved densities of states 
(DOS) for the LaNiOs structures reported in Table [TT] with 
the LSDA+£7 exchange-correlation functional. 



occupy the 2b Wyckoff position with 3 site symmetry, all 
Ni-0 bond lengths are required to be equivalent. We thus 
observe that our calculated Ni-0 bond lengths respond 
proportional to the Ni06 rotation angle (Table |Tl|). 

We next examine the change in the electronic structure 
with increasing correlation to evaluate how the electronic 
states around the Fermi level are modified. We show 
in Fig. [6] the electronic densities-of-states (DOS) as a 
function of the Hubbard U value obtained by LSDA+/7 
calculations with the VASP code. Consistent with earlier 
band structures calculations on bulk LaNiOs with the 
LSDA (U = eV) functionalJ^EUwe find a non-magnetic 
ground-state with localized Ni 3d states peaked centered 
1.0 eV below the Fermi level (Ep). A set of delocalized Ni 
3d states cross Ep, while the O 2p states are distributed 
throughout the entire valence band. 

Here, we apply a Hubbard U term on the Ni-d states 
and anticipate the on-site Coulomb interaction to localize 
the itinerant Ni-d electrons at Ep and stabilize a fer- 
romagnetic (FM) spin configuration. Indeed, we find a 
half-metallic FM ground-state (1 /ib/Lu.). Note, this con- 
figuration, however, has not been experimentally reported. 
For U = 3 eV, we find that the DOS at the Fermi level is 
substantially suppressed compared to the LSDA, with the 
Ni d states shifting to lower energy. Due to the increased 



TABLE II: Optimized internal structural parameters of R3c 
LaNiOs at the experimental cell parameter a and rhombo- 
hedral angle a, calculated with the LSDA+[7 functional and 
PAW pseudopotentials as implemented in VASP. 





Exp. 


LSDA+£7 (U in 
3 


eV) 
6 


X 


0.797 


0.796 0.800 


0.797 


d(Ni-O) (A) 


1.933 


1.932 1.935 


1.933 


ZNi-O-Ni (°) 


164.8 


165.2 163.9 


164.6 


Ni0 6 rotation (°) 


9.20 


8.97 9.77 


9.33 



TABLE III: Raman-active vibrational modes computed using 
the LSDA+£7 formalism at the experimental cell volume are 
compared to the measured values. All mode frequencies are 
given in wavenumbers (cm -1 ). The experimental values are 
taken from Refl631 





E 9 


E 9 




Eg 


E 9 


Exp. 




156 


209 


399 


451 


U = eV 


61.8 


155.0 


215.4 


372.6 


451.7 


U = 3 eV 


94.3 


164.8 


245.5 


387.3 


413.5 


U = 6 eV 


54.6 


156.7 


219.4 


316.0 


381.8 



rotational angle of the Ni06, we find a pseudo-gap opens 
in the majority spin DOS around -3.5 eV with reduced 
hybridization between the Ni 3d - O 2p states. With a 
further increase of the U value [Fig. [6jc)], we find that 
most of the majority spin Ni 3d electrons are shifted com- 
pletely to the bottom of the valence band around -5.5 eV. 
The states at Ep are now mainly O 2p character with a 
small contribution from the Ni 3d electrons. 

We now compute the Raman active lattice modes for 
the structures listed in Table [II] to explore the change of 
Raman frequencies with an increasing Hubbard U value. 
We begin by comparing the experimentally measured val- 
ues to our calculated ones (Table [m]). The best agreement 
with the experimental data is for the LSDA (U = eV) 
exchange-correlation functional. A Hubbard U value of 
3 eV overestimates the low frequency Raman mode with 
Eg symmetry (La anti-parallel displacements), while it 
underestimates the two high frequency E g modes (Ni-0 
bond bending and stretching). 

Interestingly, with the LSDA+Z7=6 eV functional 
we find an unstable (imaginary) zone-center phonon 
(246i cm -1 ) with A^ g symmetry indicating that the rhom- 
bohedral structure with the simple a~ a~ a~ tilt pattern 
is unstable. The atomic displacement pattern of the imag- 
inary mode corresponds to a three-dimensional checker- 
board arrangement of dilated and contracted octahedra — 
the so called octahedral "breathing" mode that often 
accompanies charge disproportionation reactions. 67 This 
mode, however, has not been observed in any temperature- 
dependent x-ray studies on LaNiOs. We therefore suggest 
that this correlation-induced octahedral distortion is a 
spurious artifact of using too large of an on-site Coulomb 
repulsion interaction in the LSDA+/7 calculation. We 
conclude that a Hubbard U value less than 6 eV should 
be used when simulating LaNiOs, because of the overall 
poor accuracy in the calculated Raman modes (Table 
III) despite such large Hubbard values closely reproduc- 
ing the experimental lattice parameters of LaNiOs phase 
(Table gj. 

We now explore how changes in the electron correla- 
tions modify the lattice volume by fully relaxing both the 
intern al a tom positions and the rhombohedral structure 
(Table |IV| . For the LSDA+/7 exchange-correlation poten- 
tial, we qualitatively find the same structural trends as 
in Table ITU While the LSDA underestimates the atomic 
volume, when we increase the value of U, the cell vol- 
ume increases; this can be understood as a result from 
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TABLE IV: Fully optimized structural parameters of R3c LaNiC>3 (rhombohedral setting), calculated from the various exchange- 
correlation functionals. 





Exp. 


VASP: LSDA+LT (U 


in eV) 




QE 







3 


a 



t on A 


f~\ f~\ A "DT3T7 1 


rBLsol 


X 


0.7968 


0.787 


0.795 


0.792 


0.801 


0.813 


0.807 


a (A) 


5.3837 


5.303 


5.308 


5.319 


5.324 


5.407 


5.367 


«rho (°) 


60.8596 


60.72 


60.92 


60.83 


61.39 


61.43 


61.66 


d(Ni-O) (A) 


1.933 


1.896 


1.905 


1.906 


1.923 


1.964 


1.952 


ZNi-O-Ni (°) 


164.8 


167.9 


165.6 


165.5 


163.3 


159.6 


159.6 


Ni0 6 rotation (°) 


9.2 


7.35 


8.76 


8.19 


10.08 


12.31 


11.17 


Q (A 3 /f.u.) 


56.2386 


53.57 


53.97 


54.20 


55.03 


57.71 


56.69 



enhanced electrostatic repulsion. 

We also compute the LNO ground-state structures 
with the LSDA, PBE and PBEsol functionals and norm- 
conserving pseudopotentials as implemented in QE. Typi- 
cally, the LSDA underestimates the lattice constant a, cell 
volume and Ni-0 bond length, but it closely reproduces 
the experimentally measured Ni-06 octahedral rotation 
angel 0. For the PBE functional, we find that both the 
bond length and octahedral rotation angles are overesti- 
mated. On the other hand, the PBEsol gradient-corrected 
functional corrects some of above-mentioned overestima- 
tion from PBE, although still slightly overestimating the 
rotation angle, it does provide the best agreement with 
the experimental lattice constant a, Ni-0 bond length 
and equilibrium atomic volume, among all the functionals 
we used for the structural optimization. 

We now use these ground-state structures listed in Ta- 



ble IV and calculate their Raman-active mode frequencies 
to examine the effect of lattice volume on the mode fre- 
quencies. Comparing these LSDA+/7 results (Table [V} 
to our previously calculated lattice modes obtained us- 
ing the same functionals but at the experimental volume 
(Table [m]), we find that the LSDA (U = eV) functional 
provides the best agreement with experimental data. As 
before, we find an unstable zone-centered Ni06 octahe- 
dral breathing mode in our calculations with U=6 eV. 
Interestingly, the LDSA+Hubbard U method accurately 
predicts the ground-state structural properties; however, 
it also leads to poor predictions for the lattice normal 
modes. We note that in this case, simply reproducing 
the correct ground-state atomic structure is not a suffi- 
cient criterion to evaluate the performance of a functional. 



TABLE V: Raman-active vibrational modes computed using 
the various exchange-correlation functionals at the relaxed cell 
parameters (Table IV) compared to the experimental values. 





Eg 


Eg 




Eg 


Eg 


Exp. 




156 


209 


399 


451 


U = eV 


72.5 


166.4 


196.8 


374.1 


519.4 


U = 3 eV 


73.4 


164.0 


221.4 


388.0 


399.4 


U = 6 eV 


66.0 


159.6 


201.2 


330.3 


387.3 


LSDA-PZ 


81.6 


171.0 


231.5 


398.5 


509.3 


GGA-PBE 


84.3 


168.9 


278.9 


351.1 


407.1 


PBEsol 


81.6 


169.3 


256.3 


399.7 


458.3 



This caveat is important to consider especially in the 
first-principles search for perovskite materials with large 
electron-phonon interactions. 

We also find a close relation between the predicted Ra- 
man phonon frequencies and structural parameters of the 
bulk rhombohedral LNO, especially for the A\ g and the 
two high-frequency E g modes. Typically, the A\ g mode 
frequency is sensitive to the octahedral rotation angle 
0, whereas the bond bending and breathing E g modes 
are substantially more sensitive to the predicted Ni-0 
bond length and cell volume respectively. The LSDA 
underestimates both the lattice constant a and the cell 
volume but it also overestimates the rhombohedral 
angle a r ho- This leads to a cancellation in errors and an 
overall good prediction of both and d(Ni-O) and there- 
fore calculated rotation (Ai g ) and bending (E g ) phonon 
frequencies that agree very well with the experiment. The 
LSDA functional only overestimates the stretching E g - 
mode frequency (by 10%) due to the underestimation of 
cell volume. 

Compared to the LSDA, the PBEsol functional im- 
proves predicti on o f the structural parameters, such as 
a and Q (Table |lV|) . Except for the overestimation of 0, 
PBEsol also predicts a d(Ni-O) bond lengths close to 
experiment. Therefore, as shown in Table [V| the two 
high-frequency E g modes are in good agreement with 
experiment, but the A\ g mode is overestimated by more 
than 20%. The PBE functional, in contrast, overestimates 
all of the major bulk LNO structural parameters, thus 
providing the poorest description of the Raman frequen- 
cies among the three functionals explored in the QE code. 
To summarize, both the LSDA and PBEsol functionals 
predict overall adequate Raman frequencies for LNO that 
are close to experiment; however, each functional still has 
deficiencies for specific vibrational modes, and in fact the 
LSDA values are closest to experiment overall. 

We emphasize here that the Ni06 rotation angle 0, 
the order parameter to characterize the structural phase 
transition in LaNiOa, is highly sensitive to the exchange- 
correlation functional. Even calculations with same func- 
tional, but different pseudopotentials (for example the 
LSDA calculations performed with the VASP and QE 
codes) yield values with obvious differences. Therefore 
we suggests that an accurate and comprehensive study of 
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various theoretical approximations on the descriptions of 
the octahedra rotation angles in rhombohedral perovskite 
oxides is needed. 



C. Electronic structure and experimental energy 
spectra 

In this section we compare our first-principles results 
with recent photoemission spectroscopy (PES) data to 
identify the degree of electronic correlations in rhombohe- 
dral LaNiC>3. In a single electron picture, PES measures 
the excitation energies for non-interacting electrons from 
the valence band into the continuum and therefore can 
be used as a first-order reference to single-particle DFT 
studies. We show in Figure [7] experimental PES dataP^ 
from a crystalline 20 nm LaNiC>3 film, and compare it 
with our first-principles calculated valence band DOS to 
evaluate the accuracy of our calculations in reproduc- 
ing the known electronic structure. In order to make 
a more accurate comparison, we first smear our calcu- 
lated DOS with a Gaussian function (FWHM=0.20 eV) 
to account for the experimental resolution and multiply 
by a 20 K Fermi-Dirac distribution; then we convolute 
an energy-dependent Lorentzian function [full width at 
half-max (FWHM)= Q.1\E - E F \ eV, 1 ^ where E F is the 
Fermi level] with our calculated DOS to include lifetime 
broadening effects of the photon-excited electrons. 

The experimental PES consists of four main features: 
peaks A and B are assigned to the Ni e g and ti g states, 
and the deeper c and D spectra to the O 2p dominant 
states The e g states (a) are clearly resolved as a sharp 
peak at the Fermi level and the strong t2 g peak is located 
at 1.0 eV below the Fermi edge. The O 2p states (c), 
located below the t<i g states, correspond to non-bonding 
O 2p states, as their hybridization with the Ni 3d states is 
restricted by symmetry. The O 2p bonding states (d) are 
much broader and located between -8 and -4 eV. Com- 
pared to experimental PES results, none of the exchange 
correlation functionals are able to reproduce the sharp 
spectral intensity of the e g peak at the Fermi level. Each 
method leads to an over derealization of the itinerant e g 
electrons (Fig. [7]). 

Despite the inability to reproduce peak A, the LSDA 
functional does exceedingly well in reproducing the exper- 
imental valence band features — both the correct valence 
bandwidth and energy peak positions. As shown in Fig. 
[8^a), the LSDA predicts a pronounced hybridization be- 
tween the Ni tig and O 2p states almost throughout the 
full valence band. The e g states, located at the Fermi 
level, are separated in energy from the Ni ti g orbitals. 

The LSDA+/7 functional is also able to reproduce the 
correct bandwidth; however, it fails to predict the correct 
energy peak positions. It shifts the ti g states (peak b) 
toward the bonding O 2p states (peak d) as observed 
in the orbital resolved DOS [Fig. |8jb)]. There is an ad- 
ditional shift of bonding O 2p spectral weight from the 
bottom of the valence band to around feature B which 
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FIG. 7: (Color) The experimental PES spectrum for a 20 nm 
crystalline LNO film (filled circles) from Ref 69] is compared 
to the density of states calculated with the LSDA+[7 (U=5.74 
eV), gradient-corrected (PBE) and hybrid (PBE0 and HSE) 
exchange-correlation functionals. The calculated data are 
smeared with Lorentzian and Gaussian functions and trun- 
cated with a Fermi-Dirac distribution to facilitate the compar- 
ison. See the main text for peak position assignments. 



gives the impression that the localized t^g states below 
the Fermi level are only broadened and not shifted. These 
two effects combine to give a reduced intensity in the non- 
bonding 2p states (peak c), where as experimentally they 
contribute greater intensity to the PES data. The PBE 
functional also predicts the correct delocalized electronic 
characters for LaNiOs, but underestimates the valence 
bandwidth [Fig. [8jc)]. Moreover, it predicts a ferromag- 
netic ground-state (0.56 \ib local magnetic moment per 
Ni), in contrast to the paramagnetic LSDA results. There- 
fore, the PBE functional predicts an electronic structure 
that is intermediate between the LSDA and LSDA+t/ 
methods. 

Both hybrid functionals give poor agreement between 
the calculated DOS and the experimental PES data. In 
each case the bandwidth is overestimated, with the PBE0 
(HSE) functional « 1.5 eV (1 eV) larger than the experi- 
mental results. This is in contrast to the PBE functional 
which produces a narrow (~6 eV wide) valence band [Fig. 
|8jc)] . As mentioned earlier, similar errors are also found in 
other itinerant magnetic metals when using the PBE0 and 
HSE functionals!^ The major peaks are also red-shifted 
by approximately 2 eV, i.e. PBE0 and HSE functionals 
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FIG. 8: Spin- and orbital-resolved density of states for rhombo- 
hedral LaNi0 3 , obtained from LSDA, LSDA+[7 ([7=5.74 eV), 
PBE, PBEO and HSE calculations. The Fermi energy is the 
energy zero. In each case the structure is fully relaxed accord- 
ing to the specific functional with the exception that the DOS 
for the hybrid XC-functionals are calculated with the LSDA 
ground- state atomic structure. 



shift the major Ni ti g states to between -7.5 and -5 eV 
below Ep. Similar to the LSDA+£7 approach, the hybrid 
functionals suppress the contribution of Ni states at Ep- 
This results in the removal of the Ni e g states from the 
Fermi level and shifts them to the bottom of the valence 
band. This shift of electronic states produces the unusual 
shoulder at -8 eV in the DOS calculated with the hybrid 
functionals [Fig. [T^d) and (e)]. 

The above comparison to the experimental data clearly 
indicates that the beyond-LSDA methods (the LSDA+t/ 
and hybrid density functionals) incorrectly describe the 
electronic structure of LaNi03. The origin for these dis- 
crepancies lies in the fact that the valence band is pri- 
marily composed of strongly delocalized Ni t<i g and O 2p 
states. Therefore, metallic LaNiOs is able to strongly 
screen the electron-electron interactions. This screened 
electron-electron interaction can lead to renormalization 
of the electronic states near the Fermi level, which are re- 
sponsible for the sharp e g peak observed experimentally. 70 
The accurate treatment of this phenomenon in LaNiOs 
requires a calculation using many-body DFT methods, 
i.e. quasi-particle GW, which includes the dynamically 
screened Coulomb interaction and therefore treats the 
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FIG. 9: Orbital-resolved density of states in the core energy 
region obtained from LSDA, LSDA+[7 (U =5.74 eV), PBE, 
PBEO and HSE calculations. The Fermi energy level is shifted 
to the energy zero. Note, the binding energies predicted by 
PBEO and HSE are not obtained by a simple rigid shift of the 
PBE results. 



electron screening effect in a dynamic (energy-dependent) 
waypl Our beyond-LS DA m ethods, such as LSDA+LT 
and hybrid functionals* 57 * 62 * however, are only included 
in the exchange-correlation functional, which is indepen- 
dent of the quasi-particle frequency. As such, neither the 
on-site Coulomb repulsion interaction from LDA+J7, nor 
the exact-exchange energy from hybrid density function- 
als, is able to reproduce the screening effects (and the 
above-mentioned discrepancies) present in LaNiOs. 

The inability of these beyond-LSDA methods to cap- 
ture the dynamic screening effects manifests as shifts in 
core- level binding energies. To explore this effect, we 
show in Fig. [9] the calculated orbital-resolved DOS in the 
core energy region (-35 through -10 eV) for LaNiOs. The 
discrete core level states are primarily La 5s, bp and O 2s. 
They are weakly overlapped and easily distinguishable 
by spectroscopic measurements. In Table |VI| we com- 
pare our calculated core-level binding energies to recent 
XPS data. 71 Although the LSDA method accurately re- 
produces the delocalized states in the valence region, it 
largely underestimates the binding energies of the core 
level states. In contrast, both hybrid XC-functionals sig- 
nificantly improve the description of the peak positions in 
the core region. Unlike the delocalized valence electrons, 
the core level states are weakly modulated by the screen- 
ing effect from the valence region, and therefore exhibit 
strongly localized electronic character. In this case, the 
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TABLE VI: Calculated binding energies (in eV) for the major 
energy peaks in the core region from this study are compared 
to those calculated with the PBEO (designated by the *) 
functional and experimental (Exp.) values, taken from Ref l71l 
The assignment of the peak positions follows that of Ref ]711 



Peak 


Exp. 


LSDA 


PBE 


LSDA+C/ 


HSE 


PBEO 


PBEO* 


b/c 


16.7 


14.2 


14.4 


13.6 


15.3 


16.0 


17.1 


d 


21.0 


17.6 


17.6 


17.0 


19.3 


19.8 


20.7 


e 


23.6 


18.7 


18.6 


18.0 


20.4 


21.1 


21.8 


f 


33.5 


30.5 


30.7 


30.0 


32.7 


33.5 


34.8 



hybrid fimctionals, through the addition of Fock-exchange, 
greatly reduce the self-interaction errors present in the 
LSDA and GGA functionals and therefore substantially 
improve the calculated energy spectra. We also note that 
within the LSDA+/7 framework, the Hubbard U does 
not affect the core level states since it is applied only to 
the valence Ni 3d orbitals, and therefore the core level 
states are only rigidly shifted in energy respect to those 
calculated with the LSDA functional. 

Based on our first-principles calculations, we find that 
strong hybridization between the Ni 3d and O 2p states re- 
duce the on-site <i-orbital Coulomb interaction in LaNiOa 
through enhanced screening. The conventional LSDA 
method accurately reproduces the valence band struc- 
ture and also provides the decent experimental agree- 
ment to the structural properties. In contrast, the hybrid 
exchange-correlation functionals work best for the core 
level states. We do find one subtle caveat: the Ni e g 
states experimentally exhibit dynamical correlation ef- 
fects, which we are unable to capture in either the LSDA, 
LSDA+Z7, or hybrid functionals. This enhanced spec- 
tral weight at the Fermi level has recently been repro- 
duced experimentally and through LDA+dynamical mean- 
field theory (DMFT) calculations E^HZU w here a dynamic 
(frequency-dependent) treatment of the Ni 3d electron- 
electron interactions leads to an enhanced effective mass 
and optical conductivity. 

To answer the principal question posed earlier regard- 
ing the metallic behavior in LaNiOs: We find that the 
three-dimensional topology of the perovskite structure 
with corner-connected Ni06 octahedra supports strong pd- 
hybridization and stabilizes the metallic state. Moreover, 
the relatively small octahedral rotations with rhombohe- 
dral symmetry and the weak electron-electron correlation 
effects are insufficient to disrupt the metallic state. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have studied the lattice normal 
modes and electronic properties of the correlated metal 



LaNiOs with first-principles calculations using a variety 
of exchange-correlation functionals. We examined the 
rhombohedral-to-cubic structural transition in LaNiOs 
within second-order Landau phase transition theory. We 
found that the A\ Q Raman-active mode acts as a clear 
descriptor for the magnitude of the octahedral rotations in 
rhombohedral perovskites with the a~ a~ a~ tilt pattern. 
We therefore suggest Raman spectroscopy is a plausi- 
ble route for quantifying the magnitude of octahedral 
rotations in rhombohedral perovskite oxides. 

Using a linear response method we showed that the 
correlation effects in LaNiOs originate from the Ni-e^ 
orbitals. We then proposed an orbital-dependent effec- 
tive Hubbard U value of 5.74 eV for LaNiOs to be used 
in the LSDA+L 7 " formalism. By comparing the results 
obtained from the various functionals with experimental 
spectroscopic data, we found an accurate treatment of 
the correlation effects in LaNiOs cannot be simply ob- 
tained by the LSDA-h£7 or hybrid functional methods. 
We identified that there are strong hybridization effects 
between the Ni t^g and O 2p within LaNiOs that result in 
enhanced screening capabilities and act to reduce the elec- 
tronic correlations in the e^ orbitals. Among the various 
DFT exchange-correlation functionals examined, we find 
that only the LSDA is capable of reproducing both the de- 
localized valence states and the experimentally measured 
lattice dynamics. 
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